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We investigate a dynamical mass generation mechanism for the off-diagonal gluons and ghosts in 
QO , SU(N) Yang-Mills theories, quantized in the maximal Abelian gauge. Such a mass can be seen as 

■ evidence for the Abelian dominance in that gauge, ft originates from the condensation of a mixed 

gluon-ghost operator of mass dimension two, which lowers the vacuum energy. We construct an 
y—t . effective potential for this operator by a combined use of the local composite operators technique 

^ ■ with algebraic renormalization and we discuss the gauge parameter independence of the results. We 

f"** ' also show that it is possible to connect the vacuum energy, due to the mass dimension two condensate 

£Nj I discussed here, with the non-trivial vacuum energy originating from the condensate (-Ajl), which 

CN| , has attracted much attention in the Landau gauge. 
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■ , I. INTRODUCTION. 

Ph. 

<D . 

. An unresolved problem of SU(N) Yang-Mills theory is color confinement. A physical picture that might explain 
confinement is based on the mechanism of the dual superconductivity 0, Q, according to which the low energy 
regime of QCD should be described by an effective Abelian theory in the presence of magnetic monopoles. These 
monopoles should condense, giving rise to the formation of flux tubes which confine the chromoelectric charges. 



X 



Let us provide a very short overview of the concept of Abelian gauges, which are useful in the search for 
magnetic monopoles, a crucial ingredient in the dual superconductivity picture. 

Abelian gauges. 

We recall that SU(N) has a [/(l)^ -1 subgroup, consisting of the diagonal generators. In jsj, 't Hooft proposed the 
idea of the Abelian gauges. Consider a quantity X(x), transforming in the adjoint representation of SU(N). 

X(x) -> U(x)X(x)U + {x) with U{x) e SU(N) . (1) 
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FIG. 1: A chromoelectric flux tube between a static quark-antiquark pair. 



The transformation U(x) which diagonalizcs X(x) is the one that defines the gauge. If X(x) is already diagonal, 
then clearly X(x) remains diagonal under the action of the [/(l)^ -1 subgroup. Hence, the gauge is only partially 
fixed because there is a residual Abelian gauge freedom. 

In certain space time points X4, the eigenvalues of X(x) can coincide, so that U(xi) becomes singular. These 
possible singularities give rise to the concept of (Abelian) magnetic monopoles. They have a topological meaning 
since -ki (S'C/(A r )/[/(l) Ar_1 ) ^ and we refer to 0,01 for all the necessary details. 

The dual superconductor as a mechanism behind confinement. 

Let us give a simplified picture of the dual superconductor to explain the idea. If the QCD vacuum contains 
monopoles and if these monopoles condense, there will be a dual Meissner effect which squeezes the chromoelectric 
field into a thin flux tube. This results in a linearly rising potential, V{r) — ar, between static charges, as can 
be guessed from Gauss' law, J EdS = cte or, since the main contribution is coming from the flux tube, one finds 
EAS w cte, hence V = — J Edr m cte x r 

An example of an Abelian gauge: the maximal Abelian gauge (MAG). 

Let be the Lie algebra valued connection for the gauge group SU(N), whose generators T A , satis- 
fying [T^jT 5 ] = f ABC T c ; are chosen to be antihermitean and to obey the orthonormality condition 
Tt(T a T b ) = -T F 5 AB , with A,B,C = 1, . . . , (N 2 - l). In the case of SU(N), one has T F = |. We de- 
compose the gauge field into its off-diagonal and diagonal parts, namely 

= A A T A = A a ^T a + A^T\ (2) 

where the indices i, j, ... label the N — 1 generators of the Cartan subalgebra. The remaining N(N — 1) off-diagonal 
generators will be labeled by the indices a, b, .... The field strength decomposes as 

= F A T A = F^T a + Fj H ,T i , (3) 

with the off-diagonal and diagonal parts given respectively by 

f£„ - DfAl - Df ^ + g f ahc AlAl , (4) 
F; v = d^Al-d^+gf^A^Al, 

where the covariant derivative D"' is defined with respect to the diagonal components A 1 ^ 

Df^d^ ab ~gf aU A^ . (5) 

For the Yang-Mills action one obtains 

Sym = -\j d ' x ( F ^ F " Ua + ■ ( fi ) 

The maximal Abelian gauge (MAG), introduced in corresponds to minimizing the functional 

H[A] = J d 4 x [A a ^ a ] (7) 
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One checks that TZ[A] does exhibit a residual U(1) N 1 invariance. 
The MAG can be recast into a differential form 

DfA^ = (8) 

Although we have introduced the MAG here in a functional way, it is worth mentioning that the MAG does 
correspond to the diagonalization of a certain adjoint operator, see e.g. |5|- 

The renormalizability in the continuum of the MAG was proven in la, , at the cost of introducing a quartic 
ghost interaction. The corresponding gauge fixing term turns out to be [a, 

Smag =s J d 4 x (c a (DfA b » + |6 a ) - | 5 / afc4 c a cV - ^gf abc c a c b c c ) , (9) 

where a is the MAG gauge parameter and s denotes the nilpotent BRST operator, acting as 

sA; = - (Dfc b + gf abc A\c c + gf aM A b ^) , sA^ = - (d,c l + gf * ab A«c b ) , 
sc a = gf aM c b c l + ^f abc c b c c , sc l = !lf iab c a c b , 

.sc a = b a , sc l = b l , 

sb a = , sV = . (10) 

Here c a , c 1 are the off-diagonal and the diagonal components of the Faddeev- Popov ghost field, while c a , b a are the off- 
diagonal antighost and Lagrange multiplier. We also observe that the BRST transformations I|1U|) have been obtained 
by their standard form upon projection on the off-diagonal and diagonal components of the fields. We remark that 
the MAG © can be written in the form 

Smag =83 [ d 4 x (l-A a A» a - ^c a c a ) , (11) 
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with s being the nilpotent anti-BRST transformation, acting as 

sA° = - (Dfc b + gf abc A b ^c c + gf aU A b ^) , sA^ = - (d^ + gf * ab A«c b ) , 

_ a = g fabt-b-t + Ijabc-b-c^ — * = 9 j iab-a-b ^ 

sc a = -b a + gf abc c b c c + gf abl c b c l + gf aM c b c l , sc l = -b l + gf lbc c b c c , 

sb a = -gf abc b b c c - gf aU b b t + gf abl c b b l sb l = -gf lbc b b c c . (12) 

It can be checked that s and s anticommute. 

Expression @ is easily worked out and yields 

S*mag = J d 4 x (b a (D ab A» b + |& a ) + c a D ab D^ bc c c + gc a f abl (D^A^) c l + gc a D ab (f bcd A f * c c d ) 



agf abl b a c b c l - g 2 f abl f cdl -g 1 c d A\A ilc - ^gf abc b a c b c c - jg 2 f abl f cdi c a c b c c c d 



- ^ g 2 f abc f adi c b c c c d c l - ^g 2 f abc f ade c b c c c d c e ^ . (13) 

We note that a = does in fact correspond to the "real" MAG condition, given by eq.(JHJ). However, one cannot set 
a = from the beginning since this would lead to a nonrenormalizable gauge. Some of the terms proportional to 
a would reappear due to radiative corrections, even if a = 0. See, for example, [30l |. For our purposes, this means 
that we have to keep a general throughout and leave to the end the analysis of the limit a — > 0, to recover condition JHJ. 

In order to have a complete quantization of the theory, one has to fix the residual Abelian gauge freedom by 
means of a suitable further gauge condition on the diagonal components A 1 ^ of the gauge field. A common choice for 
the Abelian gauge fixing, also adopted in the lattice papers p| || , is the Landau gauge, given by 



Sdi ag = a J d A x = J d 4 x (Vd^ 1 + c l d» (8^ + gf » 6 A£c 6 )) 



(14) 



4 



where c*,6* are the diagonal antighost and Lagrange multiplier. 
Abelian dominance. 

According to the concept of Abelian dominance, the low energy regime of QCD can be expressed solely in terms of 
Abelian degrees of freedom [jj . Lattice confirmations of the Abelian dominance can be found in 0, 0] . To our 
knowledge, there is no analytic proof of the Abelian dominance. Nevertheless, an argument that can be interpreted 
as evidence of it, is the fact that the off-diagonal gluons would attain a dynamical mass. At energies below the scale 
set by this mass, the off-diagonal gluons should decouple, and in this way one should end up with an Abelian theory 
at low energies. 

A lattice study of such an off-diagonal gluon mass reported a value of approximately 1.2GeV [f|. More re- 
cently, the off-diagonal gluon propagator was investigated numerically in reporting a similar result. 

There have been several efforts to give an analytic description of the mechanism responsible for the dynami- 
cal generation of the off-diagonal gluon mass. In [l2l llflj . a certain ghost condensate was used to construct an 
effective, off-diagonal mass. However, in 0] it was shown that the obtained mass was a tachyonic one, a fact 
confirmed later in ^^|. Another condensation, namely that of the mixed gluon-ghost operator [^Af l A fJ,a + ac a c a ) [3^ |. 
that could be responsible for the off-diagonal mass, was proposed in ^(|. That this operator should condense can be 
expected on the basis of a close analogy existing between the MAG and the renormalizable nonlinear Curci-Ferrari 
gauge 0,0]. In fact, it turns out that the mixed gluon-ghost operator can be introduced also in the Curci-Ferrari 
gauge. A detailed analysis of its condensation and of the ensuing dynamical mass generation can be found in (l9ll20|. 

Here, we shall report on the results of |2l|. It was investigated explicitly if the mass dimension two operator 
(^A a ^A tia + ac a c a ) condenses, so that a dynamical off-diagonal mass is generated in the MAG. The pathway we 
intend to follow is based on previous research in this direction in other gauges. In [2^. the local composite operator 
(LCO) technique was used to construct a renormalizable effective potential for the operator A A A^ A in the Landau 
gauge. As a consequence of (A A A^ A ) ^ 0, a dynamical mass parameter is generated |22|- The condensate (A A A^ A ) 
has attracted attention from theoretical pit , |24| as well as from the lattice side [2{| . It was shown by means of the 
algebraic renormalization technique [2(| that the LCO formalism for the condensate (A A A^ AX } is renormalizable to 
all orders of perturbation theory |2?1 ] . The same formalism was successfully employed to study the condensation of 
(^A A A^ A + ac A c A ) in the Curci-Ferrari gauge 0,H3|. We would like to note that the Landau gauge corresponds to 
a = 0. Later on, the condensation of A A A ilA was confirmed in the linear covariant gauges [28ll2^ |. which also possess 
the Landau gauge as a special case. It was proven formally that the vacuum energy does not depend on the gauge 
parameter in these gauges. As such, the linear, Curci-Ferrari and Landau gauges are all connected to each other. We 
managed to connect also the MAG with the Landau gauge, and as such with the linear and Curci-Ferrari gauges [2l| . 



II. RENORMALIZ ABILITY OF SU(N) YANG-MILLS THEORIES IN THE MAG IN THE PRESENCE 
OF THE LOCAL COMPOSITE OPERATOR (i.AJA" + ac a c a ) . 

To prove the renormalizability to all orders of perturbation theory, we shall rely on the algebraic renormalization 
formalism [2(|. In order to write down a suitable set of Ward identities, we first introduce external fields Q^", £l Ma , 
L l , L a coupled to the BRST nonlinear variations of the fields, namely 



with 



d 4 x (Dfc b + gf abc Alc c + gf ahi Al(?) - fi"* (d^ + gf iab Alc b ) 

+ L a (gf abi c b c l + ^f abc c b c c ^ +L l ^f ab c a c b ^ , (15) 



sf^ a = sQ^ = , (16) 
sL a = sV = . 



Moreover, in order to discuss the renormalizability of the gluon-ghost operator 

Omag = \A a ^ a + ac a c a , 



(17) 
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we introduce it in the starting action by means of a BRST doublet of external sources (J, A) 

sX = J , sJ = , (18) 

so that 

S LCO = sjd 4 x (x (^-A«A^ + ac a c^+C^p) (19) 

d 4 x (j Q^A"" + ac a c°^j + C y - aA6 a c a 

+ A^ Q (£>° 6 c b + gf^A^e) + aXc a (gf abl c b c i + |/ Qbc c h c c ) 

where £ is the LCO parameter accounting for the divergences present in the vacuum correlator (Omag(x)Omag{u)) , 
which are proportional to J 2 . Therefore, the complete action 

£ = 5*YM + 5MAG + Sdiag + >Scxt + SlCO > (20) 

is BRST invariant 

sE = . (21) 

As noticed in [l^.l3l|. the gluon-ghost mass operator defined in ea. l|17|) is BRST invariant on-shell. We have written 
down in [2l| all the Ward identities, which are sufficient to prove that the most general local counterterm, compatible 
with the symmetries of the model, can always be reabsorbed by means of multiplicative renormalization. As an 
interesting by-product, we have been able to establish a relation between the anomalous dimension of the gluon-ghost 
operator Omag and other, more elementary, renormalization group functions. Explicitly, it holds to all orders of 
perturbation theory that 

7 OMA o(3 2 ) = -2(^-7c.(.9 2 )), (22) 
where (3(g 2 ) = A 4 7^7 and 7c i (ff 2 ) denotes the anomalous dimension of the diagonal ghost field. 



A. The effective potential via the LCO method. 

We present here the main steps in the construction of the effective potential for a local composite operator. A 
more detailed account of the LCO formalism can be found in 133, |33| . 



To obtain the effective potential for the condensate (Omag), we set the sources f2*, f2", L a , L % and A to 
zero and consider the renormalized generating functional 



exp(-tW(J)) = / [£V]expiS(J) , 

S(J) = S YM + Smag + Sdiag + Scount + J d*x (zjJ (^Z A A a ^ a + Z a Z c ac a c a ^j + (C + 8Q (23) 

where tp denotes the relevant fields and S'count is the usual counterterm contribution, i.e. the part without the 
composite operator. The quantity 8C, is the counterterm accounting for the divergences proportional to J 2 . Using 
dimensional regularization throughout with the convention that d = 4 — e, one has the following identification 

CoJ 2 = M^(C + ^C)J 2 . (24) 
The functional W( J) obeys the renormalization group equation (RGE) 

^T^^W^^ 92 ^" 10 ^ 9 ^ I d4xJ 7J + ^9 2 >0-^W(J)=0, (25) 
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where 

v(g 2 ,0 = m^-C- (26) 

^From ea. (|24() . one finds 

77( ff 2 ,0 = 2 7OMAG ( 5 2 )C + % 2 ,a), (27) 

with 

% 2 ,«)= ^ + 2 70MAG ( 5 2 )-/3( ff 2 )A_ a7Q ( 5 2 ) ^ (5c . (28) 



Up to now, the LCO parameter £ is still an arbitrary coupling. As explained in 32, 33], simply setting £ = would 
give rise to an inhomogeneous RGE for W( J) 



^ + Pig ] W + aia{g ^ ~ 7cw {9 ] J d xJ sJ 



W{J) = 8{g\a) d 4 x—, (29) 



and a non-linear RGE for the associated effective action T for the composite operator Omag- Furthermore, mul- 
tiplicative renormalizability is lost and by varying the value of SC,, minima of the effective action can change into 
maxima or can get lost. However, £ can be made such a function of g 2 and a so that, if g 2 runs according to /3(g 2 ) 
and a according to J a (g 2 ), C(# 2 ) a ) will run according to its RGE (|2*7Jl . This is accomplished by setting £ equal to 
the solution of the differential equation 

d d \ 

/% 2 )lTT + «7a(5 2 , ol)— C(g 2 ia ) = 2 7ew ( 5 2 )C(ff 2 , a) + 5(g 2 ,a) . (30) 



dg 2 ' ' 'da 

Doing so, W(J) obeys the homogeneous renormalization group equation 

^+P{g 2 )-^ i + o lla {g 2 )^-- 10mAG {g 2 ) J d i xJ-^jW(J)=0. (31) 

To lighten the notation, we will drop the renormalization factors from now on. One will notice that there are terms 
quadratic in the source J present in W(J), obscuring the usual energy interpretation. This can be cured by removing 
the terms proportional to J 2 in the action to get a generating functional that is linear in the source, a goal easily 
achieved by inserting the following unity, 



1 = ^ J [Da] exp i J d 4 x (~ (J - Omag - £j) 



(32) 



with N the appropriate normalization factor, in ea. l|23|) to arrive at the Lagrangian 

£(A M a) = -\F; v F^ a - \f%,F»* + £ M ag + £dia g - ^ + ^S^mag - ^ (^mag) 2 

while 



exp(-iW(J)) = J [Dtp] cxp iS a (J) , (34) 



(35) 



From eas. (!23|l and 134|) . one has the following simple relation 



6W(J) 



SJ 



= - (Omag) = - ( - ) , (36) 
.7=0 \9, 



meaning that the condensate (Omag) is directly related to the expectation value of the field er, evaluated with the 
action S a — J d 4 a:£(A A1 , a). As it is obvious from ea. l|33|l . (a) ^ is sufficient to have a tree level dynamical mass for 
the off-diagonal fields. At lowest order (i.e. tree level), one finds 



TO gluon "I f' m ghost — \l a A ■ 



off-diag. _ \g a off-diag. _ / 9 a 

Co' mghost "V a Co 

Meanwhile, the diagonal degrees of freedom remain massless 
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III. GAUGE PARAMETER INDEPENDENCE OF THE VACUUM ENERGY. 



We begin this section with a few remarks on the determination of C(5 2 5 oi). From explicit calculations in perturbation 
theory, it will become clear |4Q| that the RGE functions showing up in the differential equation l|3U|) look like 

/% 2 ) - -eg 2 ~2(p g 2 + Pig 2 + ---) , 

70mag(5 2 ) = 7o(«)5 2 +7i(a)ff 4 + • • ■ , 

7a(5 2 ) = a (a)g 2 + ai(a)g 4 H , 

5(g 2 ,a) = S {a) + 5i{a)g 2 + ■ ■ ■ . 

As such, ea. <|3(J|) can be solved by expanding ((g 2 ,a) in a Laurent series in g 2 , 



(38) 



C(g 2 )a) = CoM + Cl(a) + C2(a)5 2 



(39) 



More precisely, for the first coefficients £o, £i of the expression l|39|l . one obtains 



2/3 Co + aa o~g^ = 2 7oCo + S , 
2/3i(o + aa o~^ + aai ~^ = 2 7oCi + 2 7iCo + Si 



(40) 



Notice that, in order to construct the n-loop effective potential, knowledge of the (n+l)-loop RGE functions is needed. 

The effective potential calculated with the Lagrangian (l.'i.'iH will explicitly depend on the gauge parameter a. 
The question arises concerning the vacuum energy E vaci (i.e. the effective potential evaluated at its minimum); will 
it be independent of the choice of a? Also, as it can be seen from the equations (|4U|I . each d(a) is determined 
through a first order differential equation in a. Firstly, one has to solve for Co(o0- This will introduce one arbitrary 
integration constant Co- Using the obtained value for Co(o0> one can consequently solve the first order differential 
equation for Ci( a )- This will introduce a second integration constant Ci, etc. In principle, it is possible that these 
arbitrary constants influence the vacuum energy, which would represent an unpleasant feature. Notice that the 
differential equations in a for the Q are due to the running of a in ea. H3U|) . encoded in the renormalization group 
function 7 Q (g 2 ). Assume that we would have already shown that E v&c does not depend on the choice of a. If wc then 
set a — a*, with a* a fixed point of the RGE for a at the considered order of perturbation theory, then equation 
(|3UH determining £ simplifies to 



/?(.9 2 )^C(5 2 ,«* 



270 MA o(s 2 )C(ff 2 ,a*) + % 2 ,a*) 



(41) 



since 



7« (g 2 



. 



(42) 



This will lead to simple algebraic equations for the Q(a*). Hence, no integration constants will enter the final result 
for the vacuum energy for a = a* , and since E VSbC does not depend on a, E mc will never depend on the integration 
constants, even when calculated for a general a. Hence, we can put them equal to zero from the beginning for simplicity. 



Summarizing, two questions remain. Firstly, we should prove that the value of a will not influence the ob- 
tained value for E vac . Secondly, we should show that there exists a fixed point a*. We postpone the discussion 
concerning the second question to the next section, giving a positive answer to the first one. In order to do so, let 
us reconsider the generating functional (|34|) . We have the following identification, ignoring the overall normalization 
factors 



exp(-iW(J)) = J [Dip] expiS a (J) = i J [DtpDa] expi S(J) + J 



d 4 x ( I - - Omag 



2C \9 




where S(J) and S a (J) are given respectively by ea. l|25|l . and eci. (|33|l . Obviously, 



d_l_ 

daN 



[Da] exp 



d 4 x 



-(- 

2C \g 



o 



MAG 




da 



1 = 0, 



(43) 



(44) 
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so that 



which follows directly from 



dW(J) 
da 



d xs 



+ terms oc J 



j=o 



dS{J) 
da 



ss 



d^x [ -^c a c a ) + terms oc J 



(45) 



(46) 



We see that the first term in the right hand side of (|46|l is an exact BRST variation. As such, its vacuum expectation 
value vanishes. This is the usual argument to prove the gauge parameter independence in the BRST framework 
|26|. Note that no local operator O, with sO = Omag, exists. Furthermore, extending the action of the BRST 
transformation on the cr-field by 



sa = gsOuAG = -A» a Dfc b + ab a c a - agf a 



ibi—a—b , 



-gf abc C a C b C C 

2 



one can easily check that 



if d i x£(A l ,,a) = , 



(47) 



(48) 



so that we have a BRST invariant c-action. Thus, when we consider the vacuum, corresponding to J = 0, only the 
BRST exact term in ea. <|45|) survives. The effective action T is related to W(J) through a Legendre transformation 

r = -W(J) - J d i y.J{y)^ 1 while the effective potential V{a) is defined as 



-V(a) J d 4 x = r (~ 



If c m i n is the solution of 



and hence, 



dV{a 
da 



0, then it follows from ^ ^ T 



-J, that 



O" = CTmin => J = 



(49) 



(50) 



da 



V(a) 



d 4 x 



da 



W(J) 



j=o 



or, due to ea. l|4*5|) . 



da 



= 



(51) 



(52) 



We conclude that the vacuum energy E v&c should be independent from the gauge parameter a. 

A completely analogous derivation was performed in the case of the linear gauge j^. Nevertheless, in spite 
of the previous argument, explicit results in that case showed that E vac did depend on a. In |29l ] it was argued that 
this apparent disagreement was due to a mixing of different orders of perturbation theory. We proposed a modification 
of the LCO formalism suitable circumventing this problem and obtaining a well defined gauge independent vacuum 
energy E vac , without the need of working at infinite order j2^. Instead of the action (|2"3|l . let us consider the following 
action 



S(J) — Sym + Smag + Sdiag + da 



JT(g 2 ,a)0 M AG + ^ 2 (g 2 ,a)J 2 



(53) 



where, for the moment, T{g 2 , a) is an arbitrary function of a of the form 



f(g 2 , a) = l + fo(a)g 2 + f 1 {a)g i + 0(g 6 ) 



(54) 
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and J is now the source. The generating functional becomes 

expHW(J)) = / [D4>] expiS{J) . (55) 



Taking the functional derivative of W(J) with respect to J, we obtain 

6W(J) 



S.J 

Once more, we insert unity via 



= -T{g\a){0 M A G ) ■ (56) 
5=0 



1 = 1/ [Da] exp i J ( -i _ OMAG _ frtf, a) 2 



(57) 



to arrive at the following Lagrangian 



^a) = --F* V F^ - -F; V F^ + £ MAG + £ diag - ° 2 + * ggOMAG - ^ (Omag) 2 -(58) 
^From the generating functional 

exp(-tVV(J)) = J[D<f>]expiSz(J), (59) 

S 5 (J) = f d A x(c{A^a)+J^\ . (60) 
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it follows that 

JVV(J) 



<5J 



> 5 



(?) = 5 ^( ff 2 ,a)(0 M A G ) , (61) 



The renormalizability of the action (|35l) implies that the action (|6U[) will be renormalizable too. Notice indeed that 
both actions are connected through the transformation 



J 

T{g 2 ,a) 

The tree level off-diagonal masses are now provided by 



J = • ( 62 ) 



off-diag. /S^ TO ° ff - dia s- - A a 9 — (63) 

"^iu,„, - " Co ' ghost y Co ' 



while the vacuum configuration is determined by solving the gap equation 

dV'(cr) 



do 



- , (64) 



with Via) the effective potential. Minimizing Via) will lead to a vacuum energy E vac (a) which will depend on a and 
the hitherto undetermined functions fi(a) |41|. We will determine those functions fi(a) by requiring that E mc (a) is 
a- independent. More precisely, one has 

— -^fL = => first order differential equations in a for /;(a) . (65) 
da 

Of course, in order to be able to determine the fi(a), we need an initial value for the vacuum energy -E V ac- This 
corresponds to initial conditions for the fi(a). In the case of the linear gauges, to fix the initial condition we employed 
the Landau gauge [2j|, a choice which would also be possible in case of the Curci- Ferrari gauges, since the Landau 
gauge belongs to these classes of gauges. This choice of the Landau gauge can be motivated by observing that 
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the integrated operator J d 4 xA A A fiA has a gauge invariant meaning in the Landau gauge, due to the transversality 
condition d^A^ A = 0, namely 



(VT)- 1 min / d A x 

UeSU(N) 



(A A ) U (A» A ) U 



= J d A x{A A A^ A ) in the Landau gauge , 



(66) 



with the operator on the left hand side of ea. (|66|l being gauge invariant. Moreover, the Landau gauge is also an 
all-order fixed point of the RGE for the gauge parameter in case of the linear and Curci-Ferrari gauges. At first 
glance, it could seem that it is not possible anymore to make use of the Landau gauge as initial condition in the case 
of the MAG, since the Landau gauge does not belong to the class of gauges we are currently considering. Fortunately, 
we shall be able to prove that we can use the Landau gauge as initial condition for the MAG too. This will be the 
content of the next section. 

Before turning our attention to this task, it is worth noticing that, if one would work up to infinite order, 
the expressions (|53() and l|60|l can be transformed exactly into those of (|23[) . respectively l|35|) by means of eq-lEU 
and its associated transformation 



(67) 



so that the effective potentials V{a) and V(a) are exactly the same at infinite order, and as such will give rise to the 
same, gauge parameter independent, vacuum energy. 

IV. INTERPOLATING BETWEEN THE MAG AND THE LANDAU GAUGE. 

In this section we shall introduce a generalized renormalizable gauge which interpolates between the MAG and the 
Landau gauge. This will provide a connection between these two gauges, allowing us to use the Landau gauge as 
initial condition. An example of such a generalized gauge, interpolating between the Landau and the Coulomb gauge 
was already presented in |34j . Moreover, we must realize that in the present case, we must also interpolate between 
the composite operator ^A A A^ A of the Landau gauge and the gluon-ghost operator Omag of the MAG. Although 
this seems to be a highly complicated assignment, there is an elegant way to treat it. 

Consider again the SU(N) Yang-Mills action with the MAG gauge fixing JTTJ. For the residual Abelian 
gauge freedom, we impose 

5^ ag = J d 4 x (b%A* + c l d 2 c l + gfy (gf iab A' ia c b ) + ngf lab Al (d^) c b + ng 2 f ab f cd c a c d A\A^ C 

- ngf iab A l ^ a {b b ~gP bc c c c } ) + Kgf iab A^\Dl c c c )c b + Kg 2 f aM , (68) 

where k is an additional gauge parameter. The gauge fixing l|68|l can be rewritten as a BRST exact expression 

(1 - k) s {tdpA*) + nss Q 



d 4 x 



-A],A^ 
2 ^ 



(69) 



Next, we will introduce the following generalized mass dimension two operator, 

-AlA» a + - 
2 " 2 



O = -A^A^ + -A^A"' + ac a c a , 



by means of 



LCO 



(70) 



(71) 



J d 4 x (jO + C y - a\b a c a + \A^ a D ab c b + a\c a (gf a 



KXc'duA^ + ngf \A a A^c 



with (J, A) a BRST doublet of external sources, 



sX = J , s J = 



(72) 
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As in the case of the gluon-ghost operator (|17(l . the generalized operator of ea. (|70|l turns out to be BRST invariant 
on-shell. 

Let us take a closer look at the action 

£ = Sym + S'mag + S'ding + SLco + ^ext ■ (73) 
The external source part of the action, S ex t, is the same as given in ea. (11511 . 

Also, it can be noticed that, for k — > 0, the generalized local composite operator O of eq.{7UJ) reduces to the 
composite operator Omag of the MAG, while the diagonal gauge fixing l|t)9|) reduces to the Abelian Landau gauge 
(|14fl . Said otherwise, for k — > 0, the action £' of ea. (|73|l reduces to the one we are actually interested in and which 
we have discussed in the previous sections. 

Another special case is k — * 1, a — ► 0. Then the gauge fixing terms of £' are 

Smag + Sdiag = J d^xs (-A A d»c A ) = I d*x (c A d»D AB c B + b A d»A A ) , (74) 
which is nothing else than the Landau gauge. At the same time, we also have 



1 



lim _ O = -A A A>" A , (75) 



which is the pure gluon mass operator of the Landau gauge (221 12 

From (23, we already know that the Landau gauge with the inclusion of the operator A A A^ A is renormaliz- 
able to all orders of perturbation theory. On the other hand, we have already proven the renormalizability for k = 0. 
The complete action X', as given in ea. (|73|l . is BRST invariant 

sYl = . (76) 

In [2ll |. we have written down the Ward identities of this model for n ^ and general a, and we have proven the 
renormalizability to all orders of perturbation theory. It was found that the additional gauge parameter k does not 
renormalize in an independent way, while also a generalized version of the relation l|22|l emerges 

lo( 9 2 ) = -2(^-lA9 2 )) • (77) 

Summarizing, we have constructed a renormalizable gauge that is labeled by a couple of parameters (a, n). It allows 
us to introduce a generalized composite operator O, given by ea. (|7U|) . which embodies the local operator A A A^ A of 
the Landau gauge as well as the operator Omag of the MAG. To construct the effective potential, one sets all sources 
equal to zero, except J, and introduces unity to remove the J 2 terms. A completely analogous argument as the one 
given in section III allows to conclude that the minimum value of V(a), thus E vstc , will be independent of a and n, 
essentially because the derivative with respect to a as well as with respect to k is BRST exact, up to terms in the 
source J. This independence of a and k is again only assured at infinite order in perturbation theory, so we can 
generalize the construction, proposed in section III, by making the function T of ea. (|54|l also dependent on k. The 
foregoing analysis is sufficient to make sure that we can use the Landau gauge result for E vac as the initial condition for 
the vacuum energy of the MAG. Moreover, we are now even in the position to answer the question about the existence 
of a fixed point of the RGE for the gauge parameter a, which was necessary to certify that no arbitrary constants 
would enter the results for E vac . We already mentioned that the Landau gauge, i.e. the case (a, k) = (0,1), is a 
renormalizable model |2?| . i.e. the Landau gauge is stable against radiative corrections. This can be reexpressed by 
saying that (a, k) = (0, 1) is a fixed point of the RGE for the gauge parameters, and this to all orders of perturbation 
theory 



V. NUMERICAL RESULTS FOR SU(2). 



After a quite lengthy formal construction of the LCO formalism in the case of the MAG, we are now ready to 
present explicit results. In this paper, we will restrict ourselves to the evaluation of the one-loop effective potential 
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in the case of SU(2). As renormalization scheme, we adopt the modified minimal substraction scheme (MS). Let us 
give here, for further use, the values of the one-loop anomalous dimensions of the relevant fields and couplings in the 
case of SU(2). In our conventions, one has j^EHEESl 



while 



and exploiting the relation (1221 



a result consistent with that of |3€ 



lAg 2 ) = (-3 -a) 



16tt 2 



0(g 4 ) 



7a (<? 2 ) 



Pig 2 ) 



-2a 



6\ g 2 



a J 167r 



22 g 4 
3 16tt 2 



+ 0(g 4 ) , 



0(5 6 ) 



7o M ag(s 2 ) = 



20 



- 2a 



16tt 2 



0(g 4 



(78) 
(79) 

(80) 
(81) 



The reader will notice that we have given only the 1-loop values of the anomalous dimensions, despite the 
fact that we have announced that one needs (n + l)-loop knowledge of the RGE functions to determine the n-loop 
potential. As we shall see soon, the introduction of the function ^(g 2 ^) and the use of the Landau gauge as initial 
condition allow us to determine the 1-loop results we are interested in, from the 1-loop RGE functions only. 

Let us first determine the counterterm <5£. For the generating functional W(J), we find at 1-loop ^2] 



W(J) = I d«x(-(( + 50 Y 



and employing 

In det 

with 



5 In det [S ab (d 2 + aj)\ - - In det 



s ab [ (d 2 + j) 9lat - ( i - - ) 



s ab (d 2 + j) 9liv 



S aa [(d - l)trln (d 2 + J)+ trln (d 2 + aj)} 



6 aa = N(N - 1) = 2 for N = 2 
one can calculate the divergent part of ea. (|82|l . 



W(J) = / (fx 



J 2 

- sc Y-i6, 2 



3 o 1 1 9 9 1 

3 -~ T7T^ a J 
e lbn z e 



1 2 T2 1 

87H e 



Consequently, 



1 



(a 2 -3) 1 - + 0(g 2 ) 



Next, we can compute the RGE function S(g 2 , a) from ea. (|28|l . obtaining 



5(g 2 ,a) 



a 



8ir 2 



0(g 2 



(82) 



(83) 



(84) 



(85) 



(86) 



(87) 



Having determined this, we are ready to calculate Co- The differential equation (|40|l is solved by 

Co(a) =a+(9-4a + 3a 2 )C , 

with Co an integration constant. As already explained in the previous sections, we can consistently put Co = 0. 
Here, we have written it explicitly to illustrate that, if a would coincide with the 1-loop fixed point of the RGE for 
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the gauge parameter, the part proportional to Co in ea. (|88(l would drop. Indeed, the equations 9 — 4a + 3a 2 = 
and —2a + | — ^ = 0, stemming from ea. (|79|l . are the same. Moreover, we also notice that this equation has only 
complex valued solutions. Therefore, it is even more important to have made the connection between the MAG and 
the Landau gauge by embedding them in a bigger class of gauges, since then we have the fixed point, even at all 
orders. In what follows, it is understood that £o = a. 

We now have all the ingredients to construct the 1-loop effective potential Vi(a). One obtains 

Vi{a) = %l 1 H 2/o+ ^ 9 J + ^irr^-6j-3^^r( ln c^-2j ■ ^ 

It can be checked explicitly that V\ (a) obeys the renormalization group 

jjL— Vi(a) — + terms of higher order , (90) 
by using the RGE functions (|78|I - H81|) and the fact that the anomalous dimension of a is given by 



P(g 2 ) 2 d\nT(g\ a) 
-^+7o MAG (<7 ) + ^ 



To{g ) = -7TT- + 7ew(.9 ) + A* — , (91) 



which is immediately verifiable from ea. (ffcilfl . 

We now search for the vacuum configuration by minimizing V\(a) with respect to a. We will put 7i 2 = ^ 

dV, ' 



to exclude possibly large logarithms, and find two solutions of the gap equation <=p- 



= 0, namely 



a = , (92) 

2Co 



g 2 N 

V - 



16tH 



N=2 16tt 2 (2/oCo + Ci ) + o? In a - a 2 + 1 ' (93) 



The quantity y is the relevant expansion parameter, and should be sufficiently small to have a sensible expansion. The 
value for (a) corresponding to ea. H93(l can be extracted from the 1-loop coupling constant 

g 2 (jl) = — ly- . (94) 
A)hi^ 

MS 

The first solution Ij92(l corresponds to the usual, perturbative vacuum (E vac = 0), while ea. (|93|l gives rise to a 
dynamically favoured vacuum with energy 

E«* = - i - a 2 ) (m« as ) 4 , (95) 



64tt 2 



™" g = e^A m . (96) 

From ea. (|95fl . we notice that at the 1-loop approximation, a 2 < 3 must be fulfilled in order to have E vac < 0. In 
principle, the unknown function fo(a) can be determined by solving the differential equation 

, . „ j off— diag 

«-C'vac „ _ „ / off-diag\ , a f 2 o\ I off-diag\ am gluon n 

— = # 2a (™ gluon «) + 4 (a> - 3) (m gluon *j = 

3 - a 2 ( dy dy dC, , <9y <9Ci , <9y <9/ 



\cta aCo oa aCi a/o oa J 

with initial condition E vac (a) — i?^f c ndau . However, to solve ea. (19711 knowledge of d is needed. Since we are not 
interested in fo(a) itself, but rather in the value of the vacuum energy E vac , the off-diagonal mass TO° lu ~^ lag and the 
expansion parameter y, there is a more direct way to proceed, without having to solve the eq. (!97fl . Let us first give 
the Landau gauge value for E vac in the case N = 2, which can be easily obtained from [2^. l3Sj. 

R Landau _ 9 r^A 4 

vac 128tt 2 MS ■ 
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Since the construction is such that E vac (a) = i?^ c ndau , we can equally well solve 

9 11 i A 1 /~ 9\ / nff— diaeA ^ 

~ 1287T 2 

which gives the lowest order masses 



off-diag /9 ee \ off-diag / | 9 C 6 \ 

TO gluo„ = I 2 I A MS . m ghost = V« I g 3^2 I A MS . ( 10 °) 



The result (|100fl can be used to determine y. From eq. (|^S|l one easily finds 

36 

2(3-0*) 



y= b . (101) 

y 187 + 66 In, 9 



We see thus that, for the information we are currently interested in, we do not need explicit knowledge of (1 and 
fo- We want to remark that, if £1 were known, the value for y obtained in ea. (|101|l can be used to determine /o 
from eq.(jnHJl- This is a nice feature, since the possibly difficult differential equation 1)97(1 never needs to be solved in 
this fashion. Before we come to the conclusions, let us consider the limit a — > 0, corresponding to the "real" MAG 



Df A^ b = 0. One finds 



in 



3 



off-diag 

;luon — \ 2° " / ''TR '\b 

36 



187 + 66 In § 



0.168. (102) 



The relative smallness of y means that our perturbative analysis should give qualitatively meaningful results. 



VI. DISCUSSION AND CONCLUSION. 



The aim of this paper was to give analytic evidence, as expressed by ea. (|102|) . of the dynamical mass generation for 
off-diagonal gluons in Yang-Mills theory quantized in the maximal Abelian gauge. This mass can be seen as support 
for the Abelian dominance 0, El in that gauge. This result is in qualitative agreement with the lattice version 
of the MAG, were such a mass was also reported [MB- The off-diagonal lattice gluon propagator could be fitted by 
p i+ m 'i i which is in correspondence with the tree level propagator we find. We have been able to prove the existence of 
the off-diagonal mass by investigating the condensation of a mass dimension two operator, namely (^A^A IJ,a + ac a c a ). 
It was shown how a meaningful, renormalizable effective potential for this local composite operator can be constructed. 
By evaluating this potential explicitly at 1-loop order in the case of SU (2), the formation of the condensate is favoured 
since it lowers the vacuum energy. The latter does not depend on the choice of the gauge parameter a, at least if 
one would work to infinite order in perturbation theory. We have explained in short how to overcome the problem at 
finite order and gave a way to overcome it. Moreover, we have been able to interpolate between the Landau gauge 
and the MAG by unifying them in a larger class of renormalizable gauges. This observation was used to prove that 
the vacuum energy of Yang-Mills theory in the MAG due to its mass dimension two condensate should be the same 
as the vacuum energy of Yang-Mills theory in the Landau gauge with the much explored condensate (A^A^Y It is 
worth noticing that all the gauges, where a dimension two condensate provides a dynamical gluon mass parameter, 
such as the Landau gauge [23, the Curci- Ferrari gauges [2(|, the linear gauges [2{| and the MAG, can be connected 
to each other, either directly (e.g. Landau- MAG) or via the Landau gauge (e.g. MAG and linear gauges). This also 
implies that, if (A^A^ A ) ^ in the Landau gauge, the analogous condensates in the other gauges cannot vanish 
either. 



15 



Acknowledgments. 

D. Dudal would like to thank the organizers of this conference for the kind invitation to give a talk and for the 
opportunity to have many interesting discussions with other participants. 



[1] Y. Nambu, Phys. Rev. D10 (1974) 4262; 

G. 't Hooft, High Energy Physics EPS Int. Conference, Palermo 1975, ed. A. Zichichi; 

S. Mandelstam, Phys. Rept. 23 (1976) 245. 
[2] G. t Hooft, Nucl. Phys. B 190 (1981) 455. 

[3] A. S. Kronfeld, G. Schierholz and U. J. Wiese, Nucl. Phys. B 293 (1987) 461. 

[4] A. S. Kronfeld, M. L. Laursen, G. Schierholz and U. J. Wiese, Phys. Lett. B 198 (1987) 516. 

[5] K. Amemiya and H. Suganuma, Phys. Rev. D 60 (1999) 114509. 

[6] H. Min, T. Lee and P. Y. Pac, Phys. Rev. D 32 (1985) 440. 

[7] A. R. Fazio, V. E. R. Lemes, M. S. Sarandy and S. P. Sorella, Phys. Rev. D 64 (2001) 085003. 

[8] V. G. Bornyakov, M. N. Chernodub, F. V. Gubarev, S. M. Morozov and M. I. Polikarpov, Phys. Lett. B 559 (2003) 214. 
[9] Z. F. Ezawa and A. Iwazaki, Phys. Rev. D 25 (1982) 2681. 
[10] T. Suzuki and I. Yotsuyanagi, Phys. Rev. D 42 (1990) 4257. 

[11] S. Hioki, S. Kitahara, S. Kiura, Y. Matsubara, O. Miyamura, S. Ohno and T. Suzuki, Phys. Lett. B 272 (1991) 326 

[Erratum-ibid. B 281 (1992) 416]. 
[12] M. Schaden, hep-th/9909011 

[13] K. I. Kondo and T. Shinohara, Phys. Lett. B 491 (2000) 263. 

[14] D. Dudal and H. Verschelde, J. Phys. A 36 (2003) 8507. 

[15] H. Sawayanagi, Phys. Rev. D 67 (2003) 045002. 

[16] K. I. Kondo, Phys. Lett. B 514 (2001) 335. 

[17] G. Curci and R. Ferrari, Nuovo Cim. A 32 (1976) 151. 

[18] G. Curci and R. Ferrari, Phys. Lett. B 63 (1976) 91. 

[19] D. Dudal, H. Verschelde, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro, S. P. Sorella, M. Picariello, JA. Gracey, Phys. 
Lett. B 569 (2003) 57. 

[20] D. Dudal, H. Verschelde, V. E. R. Lemes, M. S. Sarandy, S. P. Sorella and M. Picariello, Annals Phys. 308 (2003) 62. 
[21] D. Dudal, J. A. Gracey, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro, S. P. Sorella and H. Verschelde, Phys. Rev. D 70 
(2004) 114038. 

[22] H. Verschelde, K. Knecht, K. Van Acoleyen and M. Vanderkelen, Phys. Lett. B 516 (2001) 307. 
[23] F. V. Gubarev, L. Stodolsky and V. I. Zakharov, Phys. Rev. Lett. 86 (2001) 2220. 
[24] F. V. Gubarev and V. I. Zakharov, Phys. Lett. B 501 (2001). 

[25] P. Boucaud, A. Le Yaouanc, J. P. Leroy, J. Micheli, O. Pene and J. Rodriguez-Quintero, Phys. Rev. D 63 (2001) 114003. 

[26] O. Piguet and S. P. Sorella, Lect. Notes Phys. M28 (1995) 1. 

[27] D. Dudal, H. Verschelde and S. P. Sorella, Phys. Lett. B 555 (2003) 126. 

[28] D. Dudal, H. Verschelde, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro, S. P. Sorella and J. A. Gracey, Phys. Lett. B 574 
(2003) 325. 

[29] D. Dudal, H. Verschelde, J. A. Gracey, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro and S. P. Sorella, JHEP 0401 (2004) 
044. 

[30] K. I. Kondo, Phys. Rev. D 57 (1998) 7467. 

[31] K. I. Kondo, T. Murakami, T. Shinohara and T. Imai, Phys. Rev. D 65 (2002) 085034. 

[32] H. Verschelde, Phys. Lett. B 351 (1995) 242. 

[33] K. Knecht and H. Verschelde, Phys. Rev. D 64 (2001) 085006. 

[34] L. Baulieu and D. Zwanziger, Nucl. Phys. B 548 (1999) 527. 

[35] T. Shinohara, T. Imai and K. I. Kondo, Int. J. Mod. Phys. A 18 (2003) 5733. 

[36] U. Ellwanger and N. Wschebor, Int. J. Mod. Phys. A 18 (2003) 1595. 

[37] K. I. Kondo, hep-th/0303251 

[38] R. E. Browne and J. A. Gracey, JHEP 0311 (2003) 029. 

[39] The index a runs only over the N(N — 1) off-diagonal generators. 

[40] See section V. 

[41] At first order, E vax: will depend on fo(ct), at second order on fo(a) and fi(ct), etc. 
[42] We will do the transformation of W(J) to W(J) only at the end. 



